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Abstract 

In this letter we study the local persistence of the two-dimensional Blume- Capel Model by 
extension of the concept of Glauber dynamics. We verify that for any value of the ratio a = D/J 
between anisotropy D and exchange J the persistence shows a power law behavior. In particular 
for a < we find a persistence exponent Q\ = 0.2096(13), i.e. in the Ising universality class. For 
a > (a / 1) we observe the occurrence of blocking. 
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Among the many relevant quantities of interest in the modelling of nonequilibrium dy- 

n 

namics of spin systems such as first time passage which has been extensively discussed 
in the literature due to the appearance of non-trivial exponents in the power law behavior 
of the first return probability as function of time, another quantity of interest is that of 
persistence, i. e. the characterization of the time it takes for a particular spin not to change 
its state from its given t = configuration. Defining P(t) as the probability that a particular 
spin will not flip up to time t, at zero temperature temperature and for the Ising and Potts 
models this quantity was shown to behave as 

p(t)~t-"', (i) 

where the persistence exponent 9i describes the non-equilibrium relaxation of the system. 
This has been determined through coarsening simulations since one expects that the fraction 
of spins which do not change up to t represents a good estimate of P(t). As a consequence 
one may introduced a global version of the concept of persistence through the quantity P g (t) 
which represents the probability that the magnetization does not change its sign from its 
t — value, as done in Q]. Recently one of the authors explored these ideas to determine the 
associated 9 g exponent in the Blume Capel Model, in particular its behavior at the critical 
and tricritical points j^. It was found that the exponent shows an abrupt change as one 
goes from the critical points (6 g ~ 0.23, Ising universality class) to the tricritical one. 

The purpose of this letter is to extend the Glauber dynamics to the Blume Capel Model 
and to study the influence of the anisotropy on the exponent d\. For the sake of clarity, we 
start with a brief overview of the Ising model. 

At T = the dynamics can be greatly affected by local blocking configurations, and 
the energy necessary to overcome them might be so high as to render the system static. 
On the other hand, for T ^ it becomes difficult to define domains because one might 
not be able to distinguish between true domains and spin flips due to thermal fluctuations. 
Nonetheless for the Ising Model these difficulties can be overcome and a power law a power 
law decay for the fraction of persistent spins was ascertained in the whole low temperature 
phase j5|. The calculated values of the exponent were 6i = 0.22, 0.22 and 0.29 for T = 0, 
T = T c /3 and T = 2T c /3 respectively. Furthermore for T > T c an exponential cutoff was 
verified. Indeed, using a natural definition of persistence, the blocking effects is sensible 
to any temperature value, and power law only happens exactly in T = 0. The dynamics 
is implemented as follows: define the excitation energy AE associated to the transition 
0i — > — o~i as AE = —2JaiSi, with Si equal to the sum of nearest neighbors to <7j. If AE < 
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then the transition occurs with probability 1. For AE = it occurs with probability 1/2 
and the transition does not occur if AE > 0. 

To extend these ideas to the Blume-Capel Model we start out with the Hamiltonian 



where <7j G {—1, 0, 1}, J > and D represents an anisotropy. As we shall see the behavior of 
the persistence exponent is strongly dependent on the value of a = D / J yielding a nontrivial 
extension of the Ising Model. To implement the Glauber dynamics in the ground state at 
(T = 0) we consider the excitation energy (in units of J ) for the transition <jj — > <Jf as 
follows: 



with Si G [—4, 4] PI Z. To extend the T = single spin Glauber dynamics from a 2-state to a 
3— state model is not straightforward and some rules have to be introduced. Let the energy 
differences in the transition from <jj to any other two spin states aj or ah be represented by 
AE\ and AE 2 respectively. There are six possibilities to consider: 

1. If AE h2 < AE 2tl < 0, then a; -> a 1 / 2 ; 

2. If AE h2 < and AE 2 ,i > 0, then a { -> a)' 2 ; 

3. If AE\ t2 = and AE 2) x > 0, then <7j — > o^' 2 with probability = 1/2; 

4. If AE 1 = AE 2 < , Oi — > a\ or a { — > aj with probability = 1/2; 

5. If AE% = AE 2 = 0,0", goes to any of the three states with probability = 1/3; 

6. If AEi > and AE 2 > 0, the transition does not occur. 

With these rules Monte Carlo simulations were performed at T = for anixi square 
lattice with L = 160. Since the sample number N s does not play an important role because 
of the absence of significant statistical fluctuations (see e.g. |6j) we used N s = 200 with 
1000 MC steps. After some exploratory simulations a total of 60 different values of a 
within a G [—3, 3] where chosen. 

After the 300th MC step convergence to power laws are found, as can be seen in figure 
1. For a < we obtained the exponent 9\ = 0.214(3) by measuring directly the slope in the 




(2) 



AEj J = (a f - <Ti) [afa + a t ) - 



(3) 
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interval [300, 1000], with error bars obtained using 5 different bins. A more precise estimate 
can be obtained using the definition of the effective exponent (local slope): 

m = ^4^k (4) 

v ; lnr 9i{t/r) v ; 
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FIG. 1: The numerical estimate of the probability P(t) in Q as described in the text. 

For r = 10 and t = 1000 we obtain 6\ = 0.2096(13). The same power law behavior 
was observed for a = and a — 1, however with slightly different exponents, as can be 
seen in table For positive values of a 1) the system shows blocking, as can seen 

a < a = a = 1 
0i 0.2096(13) 0.1964(34) 0.1993(21) 

TABLE I: Local persistence exponent values of Blume Capel Model 

in figure 2. However a distinction must be made: in a € ]0, 1[ the persistence has a fast 
decay and reaches a constant value P(t — > oo) ~ 0.322, while for a G ]1,2] the value is 
P(t -> oo) ~ 0.246. For a = 2, P(t oo) ~ 0.287) and when a > 2 there is full blocking 
with P(t — > oo) frozen to the value ~ 0.333 (initial mean fraction of spins is null). In d = 2 
and T = this is due to the fact that for a > 2 the predominant phase is a configuration 
with all spins zero. 

Our values for a < are in agreement with more recent results for the Ising Model, 
namely 

Bi = 0.209(2) (5) 
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FIG. 2: Blocking effects for positive values of a (7^ 1). 

0. The point a = separates two distinct regions: a power law behavior for negative values 
of a and a blocking phase for positive values of this ratio. At a = we have a power law 
behavior but with a different exponent. 

Another interesting behavior is observed for a — 1, where a robust power law separates 
two different blocking "phases". The point a = 2 also divides two regions of distinct blocking 
behavior: the 1 < a < 2 region and the full-blocking region a > 2. 

To conclude, our results show that nontrivial behavior in stochastic persistence can appear 
as one extends the Ising Model to allow for higher " spins" and anisotropy effects as measured 
by our parameter a. In particular for a < one has a pure 2-state Ising behavior; the same 
is not observed for different values of a, as discussed in the text. It would be interesting 
to test these ideas for different stochastic systems in the hope that there also nontrivial 
behavior might be found. 

Acknowledgements 

S.R.D. would like to thank the Mathematical Physics Group at the University of Queens- 
land for their hospitality. 



5 



References 



[1] W. Feller, " An Introduction to Probability Theory and Its Applications", 1, John Whiley &: 
Sons (1968) 

[2] B. Dcrrida, A. J. Bray, C. Godreche, J. Phys. A 27, L357 (1994) 

[3] D. Stauffer, J. Phys. A 27, 5029 (1994) 

[4] S. N. Majumdar, A. J. Bray, S. J. Cornell, C. Sire 77 , 3704 (1996) 

[5] B. Derrida, Phys. Rev. E 55, 3705 (1997) 

[6] S. Jain, Phys. Rev. E 59 42493 (1999) 

[7] S. Jain, H. Flynn, J. Phys. A 33 8383 (2000) 

[8] R. da Silva, N. A. Alves, J. R. Drugowich de Feh'cio, Phys. Rev. E 67, 057102 (2003) 



6 



